Coherent transfer of Cooper pairs by a movable grain 
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A coherent hybrid of states with different number of Cooper pairs can be built in a superconductor 
grain as a result of periodically repeated discrete encounters with bulk superconductor leads. As 
a direct manifestation of such states a non-dissipative current depending on the phase difference 
between the leads can be measured. 
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I. INTRODUCTION 



Quantum devices allowing coherent states are in fo- 
cus of intense experimental and theoretical studies, in 
particular in connection with the problem of quantum 
computation. A promising example of such systems is a 
set of superconducting circuits connected by nanometer 
size contacts. Phase coherence between two supercon- 
ductors is maintained if electron (Cooper pair) exchange 
between them is allowed. A quantum state resulting from 
such coupling is a coherent superposition of states with 
different number of Cooper pairs in the linked supercon- 
ductors. As is well known, transfer of Cooper pairs is a 
coherent process causing a ground-state non-dissipative 
current (this is just the Josephson effect ||^). 

In recent experiments Nakamura et al. 0] demon- 
strated that a quantum hybrid of two charge states can be 
created by Josephson tunneling of Cooper pairs between 
a bulk lead and a nanometer-sized superconducting grain, 
commonly referred to as a single-Cooper-pair-box. One 
of the main results of Ref. [Q is that this hybrid state 
remains coherent for a rather long time, > 20 ns. 

Nanometer-size structures in which charge transfer 
is coupled to mechanical motion have also recently re- 
ceived considerable attention. A beautiful example of 
charge transport through a Ceo molecule accompanied 
by its quantum mechanical center-of-mass motion has 
been demonstrated experimentally [Q while another im- 
plementation using a macroscopic pendulum to mediate 
current between two spatially separated electrodes was 
provided by Tuominen et al. Q . In the latter system the 
center-of-mass motion serves as a shuttle which carries 
electrons, as predicted in Ref. [||. 

The problem we address here is whether coher- 
ent charge transport is possible by using a movable 
single-Cooper-pair-box oscillating between the two leads 
thereby providing Cooper pair exchange. The Josephson 
current associated with this new mechanism for Joseph- 
son coupling is then presented. 



II. MODEL SYSTEM 

First, to demonstrate this phenomenon let us consider 
a superconducting grain initially in a state \n) with a 
fixed number of Cooper pairs n which is brought into 
contact with a bulk lead and then removed from the lat- 
ter. Since the Josephson coupling decays strongly with 
the distance, the grain finally becomes isolated in a state. 



\n). The phases 0„ will then depend 



on the superconducting phase of the bulk superconduc- 
tor. Below we will refer to such a state as a Josephson 
Hybrid ( JH) . If a grain initially in a hybrid state passes 
by the lead it acquires a charge which depends on the 
coefficients c„ e"*'*", most importantly on the phases 
as well as on the phase on the lead. 

Now let us consider a set of leads and allow the grain 
to pass the leads sequentially. Since between the leads 
the grain is decoupled from the environment, the coeffi- 
cients c„ will remain constant while the phases 0„ evolve 
according to = fi^"^ {n\Hc\n) where He is the Hamilto- 
nian of an isolated grain. Each interaction with a lead 
causes a phase-dependent charge transfer. The resulting 
state as well as the total charge transfer is determined by 
the whole "history" of the grain's motion. 

An important issue is that the JH can be created only if 
the energy differences between states with different num- 
bers of Cooper pairs are smaller than the Josephson en- 
ergy Ej. This can be achieved even if Ej tends to zero 
by gate electrodes which induce charges on the grain. If 
the gate voltage is properly chosen then the states with 
two subsequent numbers of Cooper pairs on the grain 
become equal. In other words, the Coulomb blockade 
of tunneling is lifted, as in the case of a single-electron 
transistor. As a result, the JH consists of only two states 
with subsequent numbers of Cooper pairs. 

As the grain is removed from the lead the problem 
of how long the "memory" of a proximity between two 
superconductors persists arises. This question was first 
addressed by Leggett in Ref. who considered the mo- 
tion of a superconducting grain between superconducting 
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leads. There it was assumed that the grain motion in the 
contact regions is slow enough to reach a local thermal 
equilibrium between the grain and the nearest lead. As 
sources of destruction of such a memory he discussed 
both external random fields and the so-called internal 
mutual dephasing due to relaxation to the equilibrium 
during subsequent contacts to the leads. In contrast to 
that approach, our system is intrinsically non-adiabatic 
in the contact regions. Therefore the equilibrium is not 
reached during the contacts, the system evolves according 
to the laws of quantum mechanics, and intrinsic dynam- 
ics does not occur. 

Another important difference is that we consider the 
case of strong Coulomb blockade^ when the charging en- 
ergy Ec is much larger than the Josephson coupling en- 
ergy Ej. As a result, the hybrid state is constituted of 
only two charge states differing by one Cooper pair as 
mentioned above. In this way we avoid the decoherence 
due to quantum beats between many quantum states dis- 
cussed in Ref. m . 



that A/h. In the Eqs. @) and @, C{x) is the mutual 
capacitance, while Ej{x) is the Josephson coupling en- 
ergy. Only the states with even number of electrons on 
the grain are taken into account. As it is known ||], that 
requires the inequality A > e^/C. 




III. CALCULATIONS 

For concrete calculations let us consider charge trans- 
fer between two bulk leads due to repeated alternate con- 
tacts of the grain with those leads when the grain moves 



A similar system 
10|. If we prepare 



periodically in time as shown in Fig. 
was recently realized experimentally 
the system in a state such that after one period its state 
differs from the previous one only by a phase factor, then 
the charge transfer is periodic in time. Its magnitude as 
well as the direction are determined by the phase differ- 
ence between the leads, as well as by the phase acquired 
in course of free motion. This is a new mechanism for 
Josephson charge transfer due to Cooper pair shuttling 
which is a direct manifestation of coherent transport. 

Our model Hamiltonian as a function of the grain po- 
sition X can be expressed as a sum of the electrostatic 
part 



He = [eV2C(x)] [2n + Q{x)/e] 



(1) 



and Josephson coupling with the left (i) and right {R) 
leads, 



Hj = - J2 E}ix) cosies 



(2) 



s=L,R 



Here n is the operator of the Cooper pair number on 
the grain, Q{x) is the charge induced on the grain by 
the gate, are the order parameter phase of the s:th 
lead while $ is the order parameter phase operator of 
the grain. The operators h and $ obey the commuta- 
tion relation [n, $] = i as far as we restrict ourselves by 
the ground state. To make that possible we require the 
typical frequency of the grain oscillations to be much less 



FIG. 1. Schematic of the specific system described in the 
text. A superconducting grain executes periodic motion be- 
tween two superconducting bulk leads. The presence of the 
gate ensures that the Coulomb blockade is lifted during the 
contacts between the grain and superconductor which allows 
for the grain to be in the Josephson hybrid state. As the grain 
moves between the leads Cooper pairs are shuttled between 
them creating a DC-current through the structure. 



The energy of the Josephson coupling between the 
grain and the electrodes depends strongly on position and 
is given by Ej'^{x) = Eoexp{- \6xl,r\ /I) where Sxl,r 
is the distance between the grain and the right/left con- 
tact and I is the characteristic tunneling length. 

The dynamics of the system are governed by the evo- 
lution operator 



U{t2;ti) =fexp{-ih-^ / dt'H[xit')]) 



(3) 



For periodic grain motion x{t + T) = x{t) it is con- 
venient to consider the periodic eigenstates \a) at xa = 
x{tA) (cf. with Fig. H) 



U{tA + T-tA)\a) = e-'^- \a) 



(4) 



The Liouville equation for the statistical operator p{t) is 
chosen in the form 

dp{t)/dt^^ih-\H{t),p{t)]~mm)-Mt)] ■ (5) 

The last item in the right-hand side of Eq. (H) allows for 
the relaxation to the adiabatic equilibrium distribution. 



Z"^ exp 



[-PH[x{t)] 
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Here (3 is the inverse temperature while Z = 
Tr exp (^—l3H[x{t)]^ . We specify the relaxation as v{t) = 

i/Q exp [— |fe(i)| /Z]. One needs to introduce the relax- 
ation in order to "forget" the initial conditions to Eq. (||). 
Physically, the relaxation is caused by incoherent charge 
transfer, such as quasiparticle tunneling. In the follow- 
ing we will assume that the relaxation is infinitely slow. 
Then it does not enter the results explicitly. However, it 
determines the populations through the balance equation 
(^. The latter has the simplest form in the representa- 
tion of the evolution operator, Eq. (Q). In this case p is 
diagonal in the \a) representation 



p„„ = v{t) {a\U{T-t)p{t)W{T-t) \a)/V 



(6) 



where A{t) = T^^ j dtA{t) is the time average over the 



period T of the grain's motion. Equation can be di- 
rectly obtained as a time average of Eq. (j^) . 




FIG. 2. Schematical trajectory for the system. As the grain 
passes the leads the Hamiltonian is completely dominated by 
the respective Josephson term Hj'^ while away from the leads 
the charging energy part He dominates. This implies that the 
time evolution operator can be factored into parts as in Eq. 
(I). 

Knowing the eigenstates @ one can easily calculate 
the electric current 



I ^2efJ2 Pamela. . f^T^' 



(7) 



where is the charge carried through the cross section 
xa (cf. with Fig. ^) by the state \a) during one cycle and 
is given by the general relation 

= {a\h- U\tB; tA)fiU{tB;tA) \a) = dK/d^L ■ (8) 

The dynamics of the system under consideration is most 
simple if the electrostatic energy cost per tunneling 
Cooper pair, Ec = (2e)^/2C, is larger than the Joseph- 
son coupling energy, Ej. Then the tunneling transport 
can be blocked by Coulomb effects, and resonant tun- 
neling of Cooper pairs is possible only at very specific 
values of the induced charge Q. In particular, when 
Q = Qn = — (2n-|- l)e only two states with the number of 
excess Cooper pairs at the grain which differs by one are 



possible. This leaves us inside a two-state Hilbert space 
which can be characterized by the Cooper pair number at 
the grain, |n) and |n-f 1). If at that instant the Josephson 
couphng is "on", a coherent hybrid, Ci\n) + 026*"^ |n -I- 1), 
is formed. Such a situation takes place each time the 
grain passes the lead. Consequently, the state evolution 
can be decoupled into the "scattering" regions where the 
grain is in touch with a lead and "free motion" regions 
where the grain is decoupled from the leads and its state 
evolves as that of a non-interacting system. The grain's 
trajectory is schematically shown in Fig. |^. For such mo- 
tion the time evolution operator U has the form of a 2 x 2 
matrix which can be factored into four parts pertaining 
to the two scattering events and the evolution of the state 
between the leads i.e. 



UitA+T;tA)^U.SRU+SL. 



(9) 



Here 



U± = exp(zx±a3) , x± = J dtSEc[x{t)] (10) 

describe the "free" evolution along the upper, right mov- 
ing, and the lower, left moving, paths respectively. The 
quantity SEc{x) is the difference in electrostatic energy 
for the two charge states in the hybrid while 173 is the 
Pauli matrix. The scatterings with the leads are given 
by the matrices 



Ss = exp [ids (cti cos $s -t- <T2 sin $s 



(11) 



Here are dimensionless contact times with the leads 
s = R,L defined as 'dg = J dt Ej(t). In the following 
we will assume for simplicity that 'dn = 'ds = '&■ Since 
the the exponential rapidly decays when the grain is far 
from the lead the integral can be decoupled into the in- 
tegrals over the scattering regions and the integrals over 
the regions of the free motion. Inserting Eq. (|S|) for the 
evolution operator into the expression (0) for the elec- 
tric current, one finds the relationship between the dc 
current, / and the eigenvalues of the evolution operator, 

7 = 2e/^ p,./. = e/ {dXi/d^R) {dXi/d-d) . (12) 

i 

This expression is non-trivial because it expresses the 
current in a non-equilibrium state in an equilibrium fash- 
ion. 

The eigenvalues for our two-state system can be ex- 
pressed as exp(— iAq), where Ai ~ — A2 = A is given by 
the equation 



cos A = cos^ d cos X — sin d cos $ . 



(13) 



Here we have defined x — X- +X+ and $ — ($_r — ^l) + 
(x+ ^ X-)- The dc current for the case 1?^ -i?;? = -i? is 



2e/ 



cos "d SVC? d sin $ (cos x + cos $) 
1 — (cos^ -d cos X — sin^ ■& cos <I>)2 



(14) 
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IV. DISCUSSION AND CONCLUSIONS 

The oscillating dependence of the dc current on the 
phase difference $_r, — $l is a direct consequence of the 
fact that the coherent states are controlled by the phase 
difference $, see Eq. (|l|). In this way the Cooper pair 
shuttle discussed above provides a way to prepare and 
control a coherent superposition of quantum states in a 
two-state system which is referred to as a qubit. If there 
is no phase difference, $l — $_r, but the grain's tra- 
jectory is asymmetric, x+ 7^ X-j the current still does 
not vanish. If the grain's trajectory embeds some mag- 
netic flux created by external magnetic field with vector- 
potential A(r), an extra item (27r/$o) § ■A-(r) • dx enters 
the expression for the phase difference $ which must be 
gauge- invariant. Here $0 — nhc/e is the magnetic flux 
quantum. 




FIG. 3. The magnitude of the current 7 in Eq. (|l^) in 
units of 7o = 2e/ as a function of the phases $ and x- Re- 
gions of black correspond to no current and regions of white 
to |///o| = 0.5. The direction of the current is indicated in 
the figure by signs ±. To best see the triangular structure of 
the current the contact time has been chosen to i9 = 7r/3. Jo 
contains only the fundamental frequency of the grain's motion 
and the Cooper pair charge. 

At small contact times, <C 1, 

/ « 2e/i9^ sin$ (cosx + cos$). (15) 

The smallness of this current in this case is natural be- 
cause it is the contact with the leads which creates the 
coherent state. 

To make the system under consideration realistic, one 
has to design a setup with large decoherence times. 
Firstly, the device should be carefully screened from elec- 
tromagnetic perturbations as well as from other time- 



dependent magnetic fields. Secondly, the gate electrode, 
as well as the insulating region between the gate and the 
grain, should contain as little as possible of charged dy- 
namic defects. Time-dependent switching between the 
states of those defects would produce charge fluctuations 
coupled to the charge states of the grain. It seems very 
difficult to estimate theoretically the total dephasing. For 
the estimate of the dephasing one can take the experi- 
mental results obtained by Nakamura et al. ||l| who have 
observed coherent response of a superconducting qubit 
during the time = 10~^ — 10~* s. This is a lower bound 
for our case because the grain is coupled to the leads only 
during a small part of the period of its motion. Conse- 
quently, the contact time is of the order of ~ T \JTfd 
where d is the distance between the leads. It is this quan- 
tity that should be compared to the dephasing time, r^, 
since the rest of the period the grain is decoupled from 
the leads. One can imagine another implementation of 
the Cooper pair shuttling where the alternating couplings 
between the grain and the leads are controlled by time 
dependent potential barriers rather than by mechanical 
motion. To conclude, we have demonstrated a possibil- 
ity to create a coherent quantum superposition between 
two charge states of a superconducting grain by repeated 
alternating contacts with two (or more) superconducting 
leads. This state can be controlled by the phase difference 
between the leads and monitored by a phase-dependent 
dc current though the closed circuit including the leads. 
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